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ABSTRACT 

We present deep g, i-band DECam stellar photometry of the Hercules Milky Way satellite galaxy, 
and its surrounding field, out to a radial distance of 5.4 times the tidal radius. We have identified nine 
extended stellar substructures associated with the dwarf; preferentially distributed along the major 
axis of the galaxy. Two significant over-densities lie outside the 95% confidence band for the likely 
orbital path of the galaxy and appear to be free-floating tidal debris. We estimate the luminosity of 
the new stellar substructures, and find that approximately the same amount of stellar flux is lying in 
these extended structures as inside the main body of Hercules. We also analyse the distribution of 
candidate blue-horizontal-branch stars and find agreement with the alignment of the substructures at 
a confidence level greater than 98%. Our analysis provides a quantitative demonstration that Hercules 
is a strongly tidally disrupted system, with noticeable stellar features at least 1.9 kpc away from the 
galaxy. 

Subject headings: galaxies: dwarf (galaxies:) Local Group 


1. INTRODUCTION 

Ultra-faint Milky Way (MW) satellite galaxies are the 


2012 ). 


Mateo 

1998; Simon & Geha 

2007 

McConnachie 


pressure-dominated systems are determined from their 
velocity dispersion, which assumes that the underlying 
stellar populations are in dynamic-equilibrium. How¬ 
ever, satellite galaxies interacting with their hosts can 
undergo significant tidal stirring (Lokas et al. 2012), leav¬ 
ing kinemati c sa m ples poten tial ly co ntaminated by un¬ 
bound stars (Klimentowski et al. 2007). In the case where 
a galaxy is being tidally disrupted, the mass-to-light ra¬ 
tio can be overestimated, and th is may not be apparen t 
in the galaxy’s dispersion profile (Penarrubia et al. 2008). 
It is also possible for a galaxy that has u ndergone tidal 
disru ption to retain its spherical shape (Munoz et al. 
20081 ). 


The dynamical state of a satellite galaxy is thus a 
highly relevant question. For that reason it is impor¬ 
tant to search for signs of tidal disruption not only in 
the central region of these galaxies, but in the vicinity 
around them where tidal debris may be present. The 
MW satellites provide an excellent opportunity for this 
type of investigation, as they are close enough to be re¬ 
solved int o individual stars and can be studied in great 
detai l (see Tolstoy et al.|2009 McConnachie|2012 Jerjen 
2012| for a discussion and census of local satellites)! 

Comprehensive studies of the MW satellite galaxies 
have led to numerous investigations of their role in 
galaxy formation. In the ACDM cosmological paradigm, 
galaxies form hierarchicall y through merger and accre¬ 
tion of smaller structures (IB lumcnt hal et aj. 1984} Van 
Den B osch et al.||2005 Giocoli et ai.||2(JU9 )T" Moore et al. 
(1999) demonstrated that dark matter sub-structure oc- 


curs on galactic scales, resulting in galaxy halos appear¬ 
ing as scaled versions of galaxy clusters. Mergers and ac- 
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cretion are a feature of galaxy clusters, therefore the out¬ 
skirts of larger galaxie s should also show sign s of merger 
and accretion events (Diemand et al. 2007). A rather 


striking example in this context is Andromeda (M31), 
which possesses numerous sate llite systems a nd cop ious 
substructure in its stellar halo (McConnachie et al. 2009; 


own Milky Way in the form of stellar streams (e.g. 

New- 

berg||2002; 

Belokurov et at. 2007a} 

Grillmair 2006 
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2010; Sesar et aLp013; 
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with 


dwarf (Ibata et al.||1994), and its great tidal tails strewn 


across the sky (Majcwsk i et al.|2003| Newby e t al.|2013|. 

Con trary to the predictions of ACDM, |Kroupa et al.| 
(20051 demonstrated that the distribution of the MW 


satellite population is inconsistent with that of a 
dwarf galaxy population drawn from cosmological sub¬ 
structure. Further investigation of the satellite popula¬ 
tions of both the MW and M31 has found, in both cases, 
that the satellites have a disc-like distribution with high 
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. This has led to the a 

ternative picture 


that a significant fraction of the c urrently known MW 


satellite galaxies is of tidal origin (Kroupa et al. 


2010 


Pawl owski et al.|2 014); forming from a major collision o: 


the MW and another galaxy. This theory accounts for 
many of the observed features of the dwarf galaxy popu¬ 
lation, including the high mass-to-lig ht ratios (explaine d 
as systems driven out of equilibrium (Yang et al. 2014)). 

Learning whether or not the dwarf population is largely 
comprised of systems driven out of equilibrium, or shows 
signs more indicative of hierarchical build-up, may be a 
key factor in determining the origin of the population as 
a whole. Most of the research into the MW dwarf galaxy 
population to date has been restricted to the main stellar 
body of each dwarf. However, it is in the outer regions 
where we expect to see more subtle signs of ti dal d isrup- 
tion. The obvious example is Sagittarius (Ibata et al. 
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), however, Ursa Minor ( 
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et al. 

2003}, Ursa Major 11 ( 

Simon & Geha 2007 
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ture in their outer regions, suggesting that they are not 
dynamically pristine stellar systems. 

The need for wide-held studies of the dwarf galaxy 
population is a recurring theme. Given recent advances 
in large-mosaic CCD cameras, a logical next step is to 
choose a target and explore its surrounding area in detail. 


The Hercules dwarf spheroidal galaxy (Table 


in the Sloan Digital Sky Survey (Belokurov e 


et al. 


detected 

|2007b) 


is, in this context, particularly interestin g. It is the most 


see Coleman et al. 

2007 

2012), which has 

ed to 


MW. Coleman et al. (2007) noted that, given the large 
heliocentric distance of Hercules at 132 kpc, in order for 
it to be tidally disrupted, it must be on a highly ellipti¬ 
cal orbit about the MW; consistent with the best model 


value for the perigalactic distanc e of r p — 18.5 kpc (Aden 


et al. 

cn 

0 

0 

CM 

). 

Sand et al. 

(2009) suggested that Hercules 

may 

3e em 

bedded in a stellar stream, finding a stellar 


extension to the northwest of Hercules, and several asso¬ 
ciated stellar ov er-densities alo ng the major axis out to 


35' (1.4kpc). Martin & Jin (2010) further developed 
the idea that Hercules is part of a stream, and suggested 
that a progenitor on its orbit may have been taken close 
enough to the MW to induce disruption from a bound 
dwarf galaxy into a stellar stream. They present a vi¬ 
able scenario in which this is the case, where the orbital 
path of Hercules is aligned with its direction of elonga¬ 
tion. The Hercules stellar popul ation is old (>12Gyr) 


Kirby et al.|2008 

Aden et al. 2011 

Kabrizio et al. 2014), 

and displays 

a spread in |Fe/H| ol 

more than 0.5 dex 

(Aden et al. 

2011 

Kirby et al. 

201 

L). It has also been 

shown to have a t 

ligh level ol ct-enha] 

icement ( 

Koch et al. 


large chemical enrichment contributi on from Typell su- 


pernovae progenitors. More recently, Koch et al. (20131 
reported a significant deficiency in neutron-capture el- 
ements and suggested that the chemical enrichment of 
Hercules was governed by massive stars, coupled with a 
low star formation efficiency. All of this indicates Her¬ 
cules has a rather colourful sta r for mati on history. 


With the excepti on of Sand et al. (2009), and more re¬ 
cently Fabrizio et al. (2014p, all observations of Hercules 


have focussed on its mam stellar body (e.g. Coleman 


et al.|[2007 Aden et ah|2009a 2011). Sand et aT~( 2009| ) 


made observations along the semi-ma.ior axis, with five 


23'x 23' fields, and Fabrizio et al. (2014) made similar 
observations with a 2 x 2 arrangement of 23' x 23' fields 
centred on Hercules. In this paper, we make use of the 
large field-of-view afforded by the Dark Energy Camera 
(DECam) on the 4m Blanco telescope at Cerro Tololo 
in Chile. We present wide-held observations (see Fig¬ 
ure [I]) in order to determine the true nature of Hercules’ 
elongated shape; including how far it extends into its 
surrounds, both parallel and perpendicular to the major 
axis. In Section [2j we describe the observations and data 
reduction process, including photometry and star/galaxy 
separation. We also discuss photometric calibration and 
completeness. In section[3j we describe the process of dis- 


Table 1 

Fundamental Parameters of Hercules 


Parameter 

Value 

Ref. 11 

Morphological type 

dSph 

(1) 

RA (J2000) 

16:31:02 

(1) 

DEC (J2000) 

+12:47:30 

(1) 

l 

28.7° 

(1) 

b 

36.9° 

(1) 

D 0 

132+ 12 kpc 

(2) 

M v 

-6.6 ±0.4 

(3) 

(m - M) 0 

20.6 + 0.2 

(2) 

VQ 

45.2+ 1.1km s” 1 

(4) 

a 

3.7 + 0.9 km s -1 

(4) 

€ 

0.63 + 0.02 

e 

113.8 + 0.6° 


rh 

330 pc 

(3) 

[Fe/H] 

-2.41 ± 0.04 

(5) 

Stellar mass (dynamic) 

2.6 x 1O 6 M 0 

(6) 

Mean age 

>12 Gyr 

(6) 

Mass-to-Light ratio 

103111 

(7) 

rtidal (theoretical) 

485 pc 

(7) 


( 2 ) 

(4) 

(6) 


References: (lj Belok urov et al.l (1 2007b 1), 

j 2UU7b , ( 5) ll Vlartin et Till 1 2008 1 ), 


Coleman et al. 


Aden et al. (2U09a 


band et al. 


it, ( 5) p erjen ( 21) 121) , 

, (7) | Aden et al . (2009b) 
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Figure 1. Field-of-view of our DECam observation (showing a 
smoothed representation of our complete catalogue from Section 

a . The ellipse in the centre delineates the orientation and half 
ght radiu s of Hercule s (|McC onnachie|2012|). Outlines for the five 
fields from Sand et al. (2009) are superimposed in black, while the 
four fields ifom |Fabnzio et~aL (20141 are shown as a dashed outline. 

criminating between Hercules and MW foreground stars 
via model isochrone selection. Section [4] contains our 
analysis of the shape of Hercules and the identification 
of statistically significant over-densities, while in Section 
[5] we investigate the luminosity of these over-densities 
and perform a further test of significance using blue- 
horizontal-branch star candidates as tracers. In section 
[6j we discuss these results in relation to the tidal disrup¬ 
tion of Hercules, and summarise our findings. 

2. OBSERVATIONS AND DATA REDUCTION 

Observations were carried out over two nights on 2013 
July 12 and July 13, using DECam at the CTIO 4m 
Blanco telescope, as part of observing proposal 2013A- 
0617 (PI: D. Mackey). DECam is comprised of a hexag- 
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onally shaped mosaic of 62 2Kx4K CCDs, each with a 
pixel scale of 0727/pix, providing a total field-of-view of 
3 square degrees. 

The data set is composed of a single DECam pointing, 
encompassing more than five times the estimated tidal 


radius of Hercules (485 pc or 1216, Coleman et al. 2007 
Aden et al. 2009b), with 2 sets of 4x 900s exposures 


in g , and 1 set of 11 x 600s exposures in i. Relatively 
poorer seeing in the initial set of g-band images resulted 
in a second set of exposures being taken to obtain appro¬ 
priate photometric depth, thus giving total integration 
times of 7200s and 6600s in g and i respectively. A box 
shape dither pattern was used for exposures in both fil¬ 
ters. This was to fill in the inter-chip gaps on the focal 
plane, and to avoid stars repeatedly falling on the same 
pixels (reducing the effect of bad pixel columns and CCD 
edge pixels). 

The images were redu ced using the DECam commu¬ 
nity pipeline^] (Valdes et al.|2014). This process included 
the application of a WCS solution to each image, the 
subtraction of sky-background, and the co-addition of 
images into a single image stack for each filter. The fi¬ 
nal product is a multi-extension FITS file, containing 
the stacked image-mosaic sliced into 9 separate tiles; one 
tile per FITS extension. The combined images for the 
Hercules data set have an average seeing across the field- 
of-view of 1749 for g and 1709 for i. 

The full-frame subtraction of the sky background from 
the hexagonal DECam images resulted in the corners of 
the final stacked FITS image containing no science data, 
but having non-zero pixel values. During preliminary 
photometry, this was found to cause problems with the 
background mapping for the corner tiles, where as much 
as 50% of the tile contained non-science pixels. In or¬ 
der to overc ome this issue, Weight Watchein (Marmo & 
Bertin 2008) was used, in conjunction witn the weight 
maps produced by the community pipeline, to mask out 
the non-science pixels. By using these weight maps, the 
background map was weighted to reflect the presence of 
bright stars and the different levels of exposure created 
by dithering the CCD mosaic, producing an optimised 
background map for conducting photometry. 

2.1. Photometry 

Aperture photometry was performed on each of the 
nine tiles of t he o and i images using the program Source 


Extractoi]^] (Bertin & Arnouts 


1996). Each of the corre¬ 


sponding masked weight maps from Weight Watcher were 
used by Source Extractor to produce the background 
map before source detection and photometry were car¬ 
ried out. The extraction process was run in two passes 
for each tile. A shallow first pass determined the av¬ 
erage FWHM {F) value for point sources in the image, 
where point sources were defined to be bright, circular 
detections with a Source Extractor ‘internal flag’ of 0 in¬ 
dicating isolation from near neighbours or image edges. 
This information served as input for the second, deeper 
pass, where aperture photometry was measured within 
lxF and 2 xF. Using this method separately on each 

1 http://www.ctio.noao.edu/noao/content/Dark-Energy- 
Camera-DECam 

2 http://www.astromatic.net/software/weightwatcher 

3 http://www.astromatic.net/software/sextractor 


of the nine tiles also ensured that the aperture size was 
allowed to reflect variations in the point spread function 
(PSF) across the large field-of-view. 

Spurious and non-stellar detections were eliminated 
from the Source Extractor catalogues by discarding ob¬ 
jects whose light profile appeared significantly different 
to the point-like sources in the catalogue. Different types 
of objects were identified by examining the difference in 
instrument magnitude, m 2 —mi (where mi and m 2 corre¬ 
spond to aperture sizes of lx F and 2x F respectively). 
Stars with a well defined PSF populate a narrow range 
of negative m 2 — mi values. Background galaxies how¬ 
ever, are more extended and have a larger relative flux 
encompassed in the m^ aperture than the mi aperture. 
Objects with m 2 —mi > 0 represent cosmic rays, or other 
spurious objects which were instantly discarded. Figure[2] 
shows m 2 —mi against instrumental magnitude for a sin¬ 
gle image tile in the z-band. Note that saturation occurs 
at an instrumental magnitude of ~ —12, and therefore 
detections brighter that this were excluded from analy¬ 
sis. The stars are the main vertical distribution of points 
(centred at m 2 — mi ~ —0.68), flaring out at the faint 
end. Extended sources appear as a plume to the left of 
this population. The boundaries of the main stellar dis¬ 
tribution were modelled with an exponential function to 
create a region of acceptability. Objects that fell out¬ 
side this region were excluded from the catalogue. Our 
own star/galaxy separation method was accompanied by 
a cut based on the Source Extractor star/galaxy flag. A 
loose cut was made with objects accepted as stars if the 
flag was > 0.35. This method of combining the light 
profile and Source Extractor flag was adopted in order 
to retain more stellar objects at fainter magnitudes. The 
Source Extractor classification flag becomes less reliable 
at faint magnitudes due to the ambiguity in determining 
the shape of objects at this level; experimentation re¬ 
vealed that adopting this combined methodology gained 
at least a magnitude in depth with relatively little cost 
or contamination. 

The star/galaxy separation was performed only on the 
i band image, due to the significantly poorer seeing in 
the g band. The results of this process were applied to 
the g band during the following cross matching step. 

The last step in creating the final catalogue of stellar 
sources was to cross reference the WCS solution coordi¬ 
nates, applied during the community pipeline processing, 
between the g and i catalogues. The cross referencing 
was car r ied out using the Stilts command line package 
(Taylor 2011). Only stars detected in both filters, and 
passed through the i band star/galaxy separation, were 
included. 

Once the complete catalogue of 87,000 stellar sources 
was established, it was calibrated onto the standard 
SDSS photometric system. This was achieved by first 
matching objects from the Hercules catalog ue to stars 


from SDSS data release 10 (SDSS DR10: Ahn et al. 


2014|. The photometric zero points (Z), and colour terms 


(c), where colour is defined as gi ns t — Unst 1 were deter¬ 
mined by comparing the SDSS magnitude to the instru¬ 
mental magnitude of the 2x FWHM aperture for each 
star. Figure |3] shows this comparison for each filter (each 
point on the plot represents the mean of 10 stars, first 
sorted by colour and then grouped into tens). The fol- 
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m 2 (2xFWHM) - m.! (lxFWHM) 

Figure 2. Example of the star/galaxy separation method used, 
where objects are distinguished by their light profile. m\ corre¬ 
sponds to flux measured inside an aperture of lx F and similarly, 
m 2 an aperture of 2x F. Background galaxies are seen as a plume 
to the left of the vertical column of stellar objects. 

lowing values were determined for g and i, using a least- 
squares fit: 


Z g = 32.347 ± 0.002, c g = 0.042 ± 0.002 
Z, = 32.067 ± 0.002, c, = 0.054 ± 0.002 


The large field-of-view afforded by DECam begets a 
significant variation in Galactic extinction across the 
field. It was, therefore, deemed inappropriate to use a 
single correction value for the whole Hercules catalogue. 
Instead, the extinction cor rect ion val ues {A„ _and A,) 
from SDSS DR10, based on|Schlegel et al.| (|1998|), were 
applied to the calibrated Hercules photometry. The 
specific coefficients used were A n = 3.793 E(B — V ) 
and A, = 2.086 E(B — V ) (|Stoughton et al.| [2002]). 
The values were applied by matching each star in the 
catalogue to its nearest SDSS neighbour, and using the 
corresponding extinction value. Figure [4] illustrates the 
extinction variation 0.165 < A g < 0.280 (a variation of 
approximately 5% in apparent luminosity) ac ross the 
field, and co rresponds to the extinction maps of[Schlcgcl 


et al. (19981 in the direction of Hercules. The final, fully 
calibrated and extinction corrected catalogue was used 
as the reference for the following analysis. 


2.2. Completeness 
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Figure 3. Comparison of the Hercules catalogue with the SDSS 
photometry, used to determine the offset and colour terms for each 
filter. These values in g and i respectively are: Z g = 32.347, 
c g = 0.042, Z, = 32.067, c h = 0.054. Each point and error bar 
represents the mean and standard deviation of 10 stars, first sorted 
by colour and then grouped into tens. 
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Figure 4. Map of Galactic extinction ( A g ) used for each star in 
the Hercules catalogue. Cor rection values used ar e those from the 
SDSS photometry, based on Schlegel et al.| (fl998[) ■ 


Artificial star tests were carried out in order to esti¬ 
mate the photometric depth of the data set, and deter¬ 
mine the photometric accuracy. The IRAF addstar task 
was used to add stars of random R.A. and Dec. to each 
image tile; 40 stars each in 0.5 mag intervals between 
20.5 and 27.0 mag. The previously described photom¬ 
etry pipeline was then run on each tile to recover the 
artificial stars. This process was repeated 30 times, to 
determine the mean recovery rate for each image tile. 
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21 22 23 24 25 26 27 


Figure 5. Black points show recovered artificial stars for each 0.5 
magnitude bin averaged across the field (in g and i). Th e dashed 
lines show th e best-fitting relation of the form described by|Fleming| 
|et~al.| p995| , which was used to determine the 50% completeness 
levels. 


The completeness level of the whole image was estimated 
by taking the average over the nine tiles. The 50% com¬ 
pleteness limit was ascerta ined by fitting t he int erpola- 
tion function described by Fleming et al. 1995) to the 
distribution of recovered artificial stars (shown in Fig¬ 
ure [5]), and determined to be 25.11 in g and 24.70 in i, 
respectively. 

The photometric accuracy was determined by measur¬ 
ing the variation in recovered magnitudes for each input 
magnitude. An exponential was fit to these measure¬ 
ments and used as a function for determining the photo¬ 
metric uncertainty for individual stars in the catalogue. 
This procedure was completed separately for each pho¬ 
tometric band. 


3. FOREGROUND DISCRIMINATION 

In order to detect stellar over-densities related to the 
Hercules dwarf galaxy, it is essential to disentangle its 
stellar population from that of the MW foreground. 
Selecting stars based on their position in the colour- 
magnitude diagram (CMD) is a useful way of performing 
such discrimination, and a model isochrone provides an 
excellent aid. 

To that end, the best-fitting model isochrone and 
known distance to Hercules were used to produce a mask, 
separating MW foreground stars from Hercules stars in 
the CMD. Considering the old (> 12 Gyr), metal poor 
(~ —2.41), population of Hercules, an isochrone of age 


15 Gyr, [Fe/H]= -2.49dex and |q/Fe1=0.6 dex f r om the 
D artm outh Stellar Evolution Database (Dotter et al. 
20081 was used (see Figure [6]). Rather than applying a 


straight cut along the isochrone, a graduated mask was 
applied, comprised of weights varying with proximity to 
the isochrone. This was done to make allowances for 
photometric uncertainties, and provide a more rigorous 
method of discrimination than a straight cut selection 
box. Note that although a horizontal branch is clearly 
visible in the CMD in Figure [6j it is not included in the 
mask. These stars are used as individual tracers of the 
Hercules population later in the analysis. The tip of the 
red giant branch (RGB) is not visible due to the satura¬ 
tion point of the data. However, based on the luminosity 
function of the model isochrone, there are relatively few 
member stars populating this part of the CMD and the 
extra signal this would provide against the foreground 
would be minimal. 

In order to create the mask, the isochrone was binned 
into 0.02 mag increments in g , and a Gaussian function 
was calculated along the colour range for each bin, with 
the width determined by the colour uncertainty, and an 
amplitude of 1. This gave stars lying on the isochrone 
the highest weight. The mask is shown in the right panel 
of Figure [6] The shape of the mask reflects the change 
in photometric uncertainty with magnitude, and encom¬ 
passes the stellar population seen within the half-light 
radius of Hercules. 

Once complete, the mask was applied to all stars in the 
CMD, assigning each star a weight, w, between 0 and 
1. This immediately enabled the selection of the more 
likely Hercules stars, by employing a cutoff at any cho¬ 
sen weight. By lowering the cutoff weight, the selection 
region about the isochrone could be widened, and sim¬ 
ilarly the region could be narrowed by choosing a more 
strict cutoff. This provided a flexible means by which to 
analyse the Hercules population. 


4. ANALYSIS 

The search for over-densities was undertaken in two 
main steps: the detection and formal identification of 
potential over-densities, followed by a comparison of each 
detectio n to a co ntro l region to quantify its significance. 
Sections |4.1| and [472] describe these procedures in detail. 

4.1. Detecting Over-Densities 

In order to detect over-densities of Hercules stars in 
the DECam field, the distribution of stars in our cat¬ 
alogue was binned into a normalised 2D histogram of 
R.A. and DEC. Using 40" x 40" bins, an image of 
200x177 pixels was created, where the pixel values re¬ 
flected countbin/countsampie/areabin- To maximise the 
contrast between the Galactic foreground and Hercules 
stars, two separate histograms were created (see Figure 
[7]) using the weights described in Sectional] The first con¬ 
tained only stars close to the model isochrone (w > 0.2), 
and the second contained everything else. By subtract¬ 
ing this ‘foreground’ map from the isochrone based his¬ 
togram, artificial structures and foreground fluctuations 
were removed, leaving a well defined outline of Hercules 
and several potential over-densities (right panel in Fig¬ 
ure [7}. This residual image was then convolved with a 
Gaussian filter of varying kernel size (3, 5, and 7 pix¬ 
els; smoothing over 120", 200" and 280" respectively), in 
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Figure 6. Left: Colour-magnitude diagram of stars within the half-light radius of He rcules. Cent re: Hess diagram of the same region as 
left, with 15 Gyr model isochrone from the Dartmouth Stellar Evolution Database (|Dotter et al.||2008) overlaid. Right: Isochrone mask 
based on the model isochrone. The mask provides weights to the Hercules catalog based on proximity; stars closer to the isochrone are 
weighted more heavily. Note that the tip of the red giant branch is not visible due to the saturation point of the data. 


order to identify features on different scales. 

The stellar features detected in this final smoothed his¬ 
togram were formally identified and labelled using the 
Python package scipy.ndimag <^J This package considers 
all pixels in an image above a specified threshold, and 
defines individual ‘segments’ based on adjoining pixels. 
Each segment was given a label for identification and 
later analysis. 

We performed some experimentation in order to se¬ 
lect a suitable threshold for the analysis. We calculated 
the mean and standard deviation of the pixel values in 
the smoothed histogram, and applied threshold values 
as some multiple of this standard deviation (ay). If the 
threshold was set too low, it was found that the algorithm 
returned many detections that were subsequently deter¬ 
mined not to be significant; conversely, if the threshold 
was set too high then the algorithm found only the cen¬ 
tral body of Hercules. Thresholds of lay and 1.5oy struck 
an acceptable balance between these two extremes, and 
we performed the segmentation process twice with these 
two different thresholds in order to probe the extent of 
each detected over-density. Pixels with values > 4oy 
above the mean were considered part of the main body 
of Hercules and excluded from labelling (this region was 
slightly larger than the area contained inside the half- 
light radius). We also performed the segmentation pro¬ 
cess using a threshold of 0.5oy in order to explore the 
very low-level structure of the central over-density sur¬ 
rounding the main body of Hercules. Although we do 
not report the results at this threshold in the same de¬ 
tail as for the other two thresholds (in order to maintain 
the clarity of the paper), we return to this calculation in 
Section 5. 

Using this regime, individual features were identified 
as ‘segments’ and labelled. As the kernel increased from 
120" to 200" and 280", the smaller details disappeared 
and only the larger features remained apparent. As a 
check, this process was also performed with a 60" ker¬ 
nel, at which point the level of detection ‘noise’ became 
too large. The results from the 120" kernel are therefore 
presented here as it represents the best compromise be¬ 
tween detail and noise. Figure [8] shows the over-densities 

4 http://docs.scipy.org/doc/scipy/reference/ndimage.html 


for the two thresholds, smoothed with the 120" kernel. 
Note that the labelling based on the loy threshold was 
adopted for the 1.5oy threshold, in order to compare simi¬ 
lar over-densities (e.g. OD 13 surrounding the main body 
of Hercules at the loy threshold, decomposes into three 
segments at 1.5oy, labelled OD 13.1, OD 13.2, and OD 
13.3 respectively). 

During the segmentation process, a control region was 
also defined for use in later analysis, consisting of all 
pixels in the 2D-histogram with a value less than the de¬ 
tection threshold. This ensured that no over-density was 
used in the control sample, and the sample was evenly 
distributed across the DECam field. 

Once the identification and labelling of segments were 
complete, the stars which fell inside the segment bound¬ 
aries were extracted from the catalogue, regardless of 
their isochrone weight. This ensured that the signifi¬ 
cance of each feature was analysed indiscriminately. Sim¬ 
ilarly, all stars corresponding to the control region were 
extracted, forming a sample for significance testing. 

4.2. Significance of Detections 

Once the stars of each over-density had been extracted, 
the significance of each feature was tested relative to the 
model isochrone. As well as this, a brief test was con¬ 
ducted to check if extinction or variations in complete¬ 
ness across the field could be responsible for the over¬ 
densities detected. The location of the over- dens ities was 
compared to the extinction map from section[2T]to check 
for correlations. It was found that the peak-to-peak am¬ 
plitude of the extinction had to be increased by a factor 
of three before there were even hints of a correlation. As 
this equates to approximately five times the variation in 
the derived 50% completeness level across the field, it 
was deemed that these factors do not play a significant 
role in our detection of over-densities. 

A feature was determined to be significant if a notable 
number of its stars are close to the isochrone in the CMD, 
i.e. w > 0.8. The number of stars with w > 0.8 (N w > o.s) 
was counted and recorded for each segment. The value 
of N w > o.s for each detection was t hen tested against the 
control sample defined in Section |4.1| In order to carry 
out this test, a sub-sample of stars from the control sam¬ 
ple was drawn randomly such that the number of stars 
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Figure 7. Demonstration of the binning and smoothing process. Left: 2D histogram of stars in the DECam field consistent with Hercules 
w > 0.2. Centre: the complement of the left histogram (w < 0.2) forming a ‘foreground map’. Right: difference of the two histograms, 
smoothed with a Gaussian filter of 3 pixel (120 // ) kernel size. 



Figure 8. Labelling of segments from the segmentation process 
for the 3 pixel (120 // ) kernel smoothed image. 1 at threshold de¬ 
tections are in grey, and 1.5<7t threshold detections are in black. 
The labelling scheme is such that the label of the 1.5 ct t threshold 
detections take on the label of their 1 at threshold counter-parts. 
Where there is more than one corresponding segment, labelling 
takes a decimal form (e.g. Segments 13.1, 13.2, and 13.3 are the 
1.5 <Jt threshold counter-parts to segment 13). Note that the very 
centre of Hercules is excluded from segmentation. 


in the sub-sample was equal to the number of stars in a 
given segment. N w >o,g was counted and recorded for this 
control sub-sample, similarly to the segment. This pro¬ 
cess was repeated 10,000 times. The resulting frequency 
distribution was then compared to the value recorded for 
the segment. Figure [9] illustrates the results for two of 
the over-densities detected (OD 13 and OD 22). The fre¬ 
quency distribution of the control sampling is displayed 
as a histogram, with a box and whisker plot showing the 
range of possible values the control samples could take. 
The value of N w > o.s recorded for the segment is marked 
with a black dashed line. The result for OD 13 (top) 
represents a significant detection, where N w > o.s for the 
over-density is greater than the range of values found 


for the control samples. The result for OD 22 (bottom), 
on the other hand, suggests that this over-density could 
be explained as a random fluctuation in the foreground 
population, with N w >o.8 for the over-density being quite 
close to the mean control sample value. 

A ‘significance value’ (£) was determined for each test 
result, as a means of identifying the most significant de¬ 
tections, and comparing these results to the null hypoth¬ 
esis. The value of £ was generated by calculating how far 
N w >o.8 was away from the mean of the control sample 
((N w > 0 .s ( CS ))). This was achieved by fitting a Gaus¬ 
sian to the frequency distribution of N w >o.$ values for 
the control samples (see Figure |9b, and expressing the 
distance in units of the standard (deviation (cr): 


c = 


N w >o.s(OD) - (N w > 0 . 8 (CS)) 


A detection was considered significant if £ > 2.0. The 
results of this analysis are summarised in Tables [2] and 
[3j and presented graphically in Figure [TO] Diamonds 
and circles in this figure represent detections based on 
the lay and 1.5cq thresholds respectively. There are nine 
over-densities that are significant. Several of these are 
sufficiently large that they appear significant in the larger 
5 pixel and 7 pixel convolutions. These will be discussed 
in more detail later. 

In order to test the null hypo thesi s, the significance 
analysis, as described in section |4.2[ was repeated for 
300 randomly selected segments from the control region, 
varying in size between 4-18 sqr arcmin. Figure |11| 
shows the resulting frequency distribution of £ values for 
the 300 segments. The mean and standard deviation are 
computed as -0.30 and 0.89, respectively, further sup¬ 
porting the significance of detections with £ > 2.0. 

The significance test was a lso performed on the 
‘nuggets’ reported by Sand et al. (2009) (see their Figure 
9). The majority of these nuggets appear quite small, 
varying from less than an arcminute to a few arcminutes 
in size. Using the coordinates of the nuggets from this 
figure as centre, two different sized boxes of 3' and 6' 
width were defined for each nugget in order to examine 
their extent, and ensure a large enough sample of stars 
in the nuggets for our statistical analysis. The segmenta- 
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Figure 9. Results of the control sampling process for two over¬ 
density detections (Top: OD 13, Bottom: OD 22). The histograms 
represent the frequency distribution of N w >q g values from the 
10,000 control sub-samples taken for each over-density. The black 
dashed lines represents N w >q g for each over-density. The solid 
grey lines show a Gaussian fit to each frequency distribution. The 
box plot represents the median and 25th/75th percentiles of the 
distribution, with the whiskers showing the most extreme data 
points not considered outliers, and the filled circles representing 
the outliers. Note that the outliers are shown at 50% opacity so 
underlying points may be seen. 


tion and significance testing was performed for each box 
size, mimicking the process used for detecting our over¬ 
densities. Th e re sults are summarised in Tables |d] and 
[5j and Figure [l2j Only one of these nuggets coincides 
with an over-density in our list: NGT 9 (NGT 9.2 at 3' 
width) showed £ = 2.21, while our corresponding over¬ 
density (OD 13.2) showed £ = 2.11. All other nugget 
segments have significance values considerably below our 
£ = 2 limit. Since their data was somewhat deeper than 
ours, it is unsurprising that we may not detect all of the 
nuggets they did. It should be noted that with the ex¬ 
ception of OD 13, there is little overl ap between o ur over- 
densities and the fields observed by Sand et al. (2009). 
Several of our less significant detections (£ < 2.0) have 
very small amounts of overlap with these fields. Since 


Table 2 

Significance values for the 1 at detected Hercules over-densities, 
where segment labelling corresponds to the grey outlines in 
Figure [8] 


Segment 

AT* 

N w >Q.&(OD) 

(Nu,>oACS)) 

c 

OUT 

18 


r~ 

1.93 

OD 2 

136 

17 

12 

1.52 

OD 3 

65 

10 

6 

1.75 

OD 4 

74 

10 

7 

1.42 

OD 5 

16 

2 

1 

0.73 

OD 6 

309 

44 

28 

3.18 

OD 7 

21 

1 

2 

0.13 

OD 8 

377 

48 

34 

2.50 

OD 9 

530 

64 

48 

2.45 

OD 10 

55 

6 

5 

0.62 

OD 11 

19 

3 

2 

1.15 

OD 12 

89 

9 

8 

0.59 

OD 13 

2468 

316 

224 

6.47 

OD 14 

39 

8 

3 

2.39 

OD 15 

4 

0 

0 

0.35 

OD 16 

908 

123 

82 

4.75 

OD 17 

20 

3 

2 

1.10 

OD 18 

53 

8 

5 

1.60 

OD 19 

5 

0 

0 

0.14 

OD 20 

257 

36 

23 

2.78 

OD 21 

53 

8 

5 

1.56 

OD 22 

150 

16 

13 

0.74 

OD 23 

226 

35 

20 

3.41 

OD 24 

21 

5 

2 

2.30 


Table 3 

Significance values for the 1.5<rt detected Hercules over-densities, 
where segment labelling corresponds to the black outlines in 
Figure [8] 


Segment 

Ah 

N w >o.8 {OD) 

(N w >0.8(CS)) 

c 

OD“6 

2U7~ 

34 

19 

3.74 

OD 8 

183 

27 

17 

2.73 

OD 9 

282 

37 

26 

2.34 

OD 13.1 

1019 

144 

93 

5.51 

OD 13.2 

255 

33 

23 

2.11 

OD 13.3 

16 

2 

1 

0.74 

OD 16.1 

104 

18 

9 

3.00 

OD 16.2 

168 

20 

15 

1.29 

OD 20 

170 

24 

16 

2.26 

OD 23 

144 

24 

13 

3.15 

OD 24 

13 

2 

1 

0.99 


these detections are not strong, it is unlikely that partial 
coverage would yield a detection. There is however, par¬ 
tial over-lap with OD 9 and OD 16, which we find to be 
significant. These over-densities appear less significant 
at the 1.5tJf threshold, indicating they are extremely dif¬ 
fuse. Partial coverage of such diffuse over-densities may 
cause them to go unnoticed, even in a deeper data set. 
With regard to the main halo-likc body of OD 13, there 
is remarkable cor relati on with the contours in Figure 9 
from Sand et al. (20091 (upper central panel). 


5. RESULTS 

Further to the investigation of the spatial distribution 
of the over-densities detected, we tried to estimate the lu¬ 
minosity of these new over-densities. For that purpose, 
the over-densities were subdivided into four regions. The 
first region was Hercules itself, as defined during the 
segmentation process. It covers about 1.75 times the 
area of the half-light ellipse and was excluded from our 
analysis. The second was the over-density surrounding 
Hercules, OD 13. The third was all other over-densities 
with > 3.0, and the fourth region included those over- 
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Figure 10. The significance value, C for each of the over-densities detected (correspondin g t o Tables [ 2 ] and [ 3 J. Grey diamonds denote a 
lcrt threshold detection, and black circles denote a 1.5ert detection (corresponding to Figure [8). 
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Figure 11. Histogram of possible values of £ for 300 blind seg¬ 
ments taken from regions of the Hercules field devoid of over¬ 
density detections. 


Table 4 

Significance values for ‘nuggets’ depicted by|Sand et al.||2009|, 
using a 6' X 6' box for segmentati on. The segments are shown in 
Figure |13| 


Segment 

N * 

N W>0 . S (NGT) {N W>0 . 8 (C'S)) 

c 

NOT 1 

32T“ 

32 

2y 

0.57 

NGT 2 

1022 

84 

93 

-0.92 

NGT 3 

248 

28 

22 

1.25 

NGT 4 

570 

53 

52 

0.21 

NGT 5 

248 

23 

22 

0.14 

NGT 6 

651 

52 

59 

-0.93 

NGT 7 

260 

23 

23 

-0.07 

NGT 8 

278 

30 

25 

1.03 

NGT 9 

717 

82 

65 

2.21 

NGT 10 

308 

16 

28 

-2.31 

NGT 11 

291 

27 

26 

0.18 


densities with 2.0 < ( < 3.0. For each region, the com¬ 
bined flux of all stars with w > 0.8 was calculated. An 
estimate of the foreground flux was made by summing 
the flux from stars with w > 0.8 for fields drawn ran¬ 
domly from the control region, equal in area to each of 
the over-density regions. This was performed 1000 times 
for each of the four over-density regions, in order to build 


Table 5 

Significance values for ‘nuggets’ depicted by|Sand et al. (2009), 
using a 3 1 X 3' box for segmentat ion. The segments are shown in 
Figure |13| 


Segment 

N, 

N W>0 . S (NCJT) 

(N w>0 .8 (CS)) 

c 

NGT 1 

94 

ir~ 

s~ 

0.98 

NGT 2.1 

353 

32 

32 

0.03 

NGT 2.2 

68 

2 

6 

-1.32 

NGT 3 

87 

8 

8 

0.16 

NGT 4.1 

74 

8 

7 

0.69 

NGT 4.2 

107 

12 

10 

0.83 

NGT 5 

84 

6 

7 

-0.31 

NGT 6.1 

84 

6 

7 

-0.30 

NGT 6.2 

164 

11 

15 

-0.98 

NGT 7 

101 

8 

9 

-0.33 

NGT 8 

84 

10 

7 

1.02 

NGT 9.1 

189 

20 

17 

0.78 

NGT 9.2 

95 

13 

8 

1.64 

NGT 10 

106 

7 

10 

-0.84 

NGT 11 

95 

8 

8 

-0.12 



Figure 12. Similar to Figure [To| this pl ot shows the res ults for 
segmentation based on the nuggets from Sand et al.| (|2009|). Black 
points represent the 3/ segments, and grey points represent the 6' 
segments. 

up a frequency distribution for which a statistically ro¬ 
bust mean and standard deviation could be determined. 
The mean value for each was taken to represent the fore- 
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ground over that area of sky, and subtracted from the 
flux counted for the corresponding over-density region. 
The remaining flux was considered to represent an up¬ 
per limit for the luminosity of the over-density regions. 
The extended region 2, surrounding the main body of 
Hercules (region 1), was calculated to have 52 ± 12% of 
the flux of region 1, spread over approximately 2.7 times 
the area. Similarly, Region 3 contains 42 ± 9% as much 
flux as region 1, spread over approximately 1.4 times the 
area. Region 4 contains 25 ± 8% as much flux as region 
1 spread over approximately 1.2 times the area. In total, 
the combined flux found in all significant over-densities 
(regions 2, 3, and 4) is of the same order (119 ± 29%) as 
the flux of Hercules itself, thus a significant portion of 
Hercules’ stellar population has been tidally disrupted. 

The summary of our analysis is illustrated in Figure 
[13] The new over-densities are colour-coded according 
to their significance, with the 1.5cr t detections laid over 
the 1 a t detections. We have also over-laid the outline of 
those 0.5(7t detections with £ > 2.0, in order to provide 
the broadest illust ratio n of Hercules. The segmentation 
of the nuggets from Sand et al. (2009) is shown as an over¬ 
lay of black outlined boxes. Stars consistent with the blue 
horizontal branch (BHB) of Hercules (21.0 < g 0 < 21.2 
—0.40 < (g — i )o < —0.15), are shown as open black 
stars. The elon gation of Herc ules a nd the postulated or¬ 
bital path from Martin & Jin| (20101 are demonstrated by 
a grey ellipse and black solid line respectively. We note 
that the orbital path presented by Martin & Jin (2010) 
represents a viable scenario in which Hercules has had a 
close encounter with the MW, and the orbital path was 
assumed to be aligned with the direction of elongation 
of Hercules. This is consistent with the idea that Her¬ 
cules has very recently passed peri-galacticon, not yet 
havin g flip ped its_ tid al tails in an orientation toward the 
MW (Klimentowski et al. 2009). Small grey-filled stars 
represent stars identified as potential Hercules members, 
by means of re l ative proper motions, from Table 2 of 
Fabrizio et al. (2014). This table includes all 528 of 
the initial proper motion candidates identified by |Fab- 


rizio et al. (2014), which they state contains a level of 
remaining held contamination. In order to reduce the 
number of contaminants in the sample, we preferentially 
selected stars sitting closest to the red-giant and hori¬ 
zontal branches in the CMD. The distribution of these 
239 proper motion, CMD, selected stars shows a remark¬ 
able correlation with the outline of the extended body of 
Hercules. 


Connachie 2012). An orientation angle of 113.8±0.6° was 
found, which is slightly steeper th an the published value 
of 102° ± 4° jMcCon nachle 2012). In the scenario from 
McConnachie"i2012|), where the orbital path is aligned 
with the position angle of the major axis of Hercules, 
this would indicate a steeper orbital path as plotted with 
the dash-dot line in Figure [l3| This path lies closer to 
several of the leading and trailing over-densities. To test 
the correlation between the properties of the main stellar 
body of Hercules and the distribution of the newly iden¬ 
tified stellar over-densities, we fit a line to the positions 
of all of the over-densities, using the ( values as weights 
for the fit. This line is shown in Figure [Ts] as the dash- 
dot-dot line. Confidence bands for this nt at the 68% 
and 95% level are shown by the dashed and d otted grey 
line s respe ctively. Both the orbital path from |Martin fe| 
Jin (2010), and our revised orientation of the major axis 
of Hercules, fall inside the 68% confidence band. 

A clear trend is also detected in the distribution of the 
BHB star candidates, following the elongation of Her¬ 
cules. Since these are BHB star candidates (as photo¬ 
metrically defined by the black box in the CMD in Fig¬ 
ure |6b , the expectation is that they belong to Hercules, 
and therefore provide a good tracer of its stellar popu¬ 
lation. Tests were performed to see if the distribution 
of BHB stars could be reproduced with a random selec¬ 
tion of stars from across the field. There are 21 BHB star 
candidates, however, the seven stars within the half-light 
radius were excluded from testing due to the high prob¬ 
ability of membership. Thus, 14 stars remained in the 
outer regions. Of these, six appear to coincide with (or 
sit on top of) the detected over-densities; in addition, six 
are observed to sit within the 68% confidence band and 
11 fall within the 95% band. It was found that, draw¬ 
ing a random sample of 14 stars from anywhere in the 
field, six or more coincided with an over-density 1364 
times out of 10 7 trials (0.01%). A similar test was per¬ 
formed to assess the likelihood of stars falling inside the 
68% and 95% confidence bands. It was found that six or 
more out of 14 stars he inside the 68% interval 0.09% of 
the time (out of 10 7 trials). Similarly, 11 or more out of 
14 stars fall inside the 95% confidence band 2% of the 
time (out of 10' trials). This shows that the distribution 
of BHB candidates we detect beyond Hercules’ radius is 
in agreement with the alignment of the substructures at 
a confidence level greater than 98%. 

6 . SUMMARY 


The overall distribution of the over-densities detected 

The tidal disruption of the Hercules dwarf galaxy has 

appears roughly aligned with the major axis of Hercules. 

been alluded to for some time ( 

Coleman et al.|2007 

Aden 

However, two of the significant detections (OD 6 and 

et al. 2009b Sand et al. 2009 

Martin & J in| 2010), how- 

OD 23) lie at notably different orientation angles to the 

ever the full extent to which 

it occurs was not easy to 


major axis. These over-densities could be free-floating 
tidal debris following Hercules on slightly different or¬ 
bits. As a means of furt her investigation, the astroMlj^] 


(Vanderplas et al. 2012) machine learning package was 
employed to determine the ellipticity and angular ori¬ 
entation of Hercules in our data set. Taking the cen¬ 
tral 9'x9'(containing the half-light radius of Hercules), a 
maximum likelihood bivariate normal fit was performed. 
An ellipticity of e = 0.63 ± 0.02 was found, which is in 


agreement with the published value of 0.68 ± 0.08 (Mc- 


5 http://www.astroml.org/ 


quantify due to technical limitations. Our deep, photo¬ 
metric analysis of 3 sqr deg around Hercules, with the 
Dark Energy Camera, reveals that this dwarf galaxy has 
undergone significant tidal disruption, with nine statis¬ 
tically significant over-densities extending to more than 
2 kpc from its centre. Most over-densities are distributed 
along the major axis, with all but two falling inside the 
95% confidence band of the likely orbital path. We es¬ 
timated the combined luminosity of these over-densities 
relative to the Hercules galaxy itself. The total number 
of stars in these over-densities is of the same order as 
the total number of stars in the main body of Hercules, 
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Figure 13. Summary figure displaying compiled results from this analysis. The significance of lat and I .oat detections are colour coded 
according to (1.5 at detections appear ‘inside’ the 1 oy dete ctions). Col oured contours reflect those over-densities detected at the 0.5ot 
threshold with ( > 2.0. Our segmentation of the nuggets from|Sand_et at (2009|) is shown as an overlay of black outlined boxes. Candidate 
blue horizontal branch st ars, selected from the CMD (see Figure HTp . are shown as open black stars. The elongation of Hercules and the 
orbital path proposed bylMartin & Jin (2010) are demonstrated Dy a black solid line and grey ellipse. Small grey-f illed stars represent 
stars identified as potential Hercules members, through relative proper motions, from Table 2 of|Fabrizio et al.|j2014;). The dash-dot line 
represents the inclination of Hercules in this data set based on a maximum likelihood bivariate normal fit to stars inside the central 9'x9h 
The dash-dot-dot line represents a weighted least squares fit to the position of the over-densities, using ( as weights. Confidence intervals 
for this fit at the 68 and 95 percentile limits are shown by the dashed and dotted grey lines respectively. 


highlighting the severe tidal disruption this dwarf galaxy 
has experienced in the tidal field of the Milky Way. 

Although we find littl e corre lation with the nuggets 
reported by Sand et al. (2009) (with the exception of 
the shape of the extended region surrounding the main 
body of Hercules), we find excellent correlation with the 
proper-moti on identified stars from Table 2 of |Fabrizio| 
et al.| (|2014|). Given the distinctly different method used 
by" Fabrizio et al. (20141, the correlation of these results 
with our own provides extra confidence in the position 
of our over-densities, and encourages the use of relative 
proper motion as a means of identifying membership. 

Our results provide useful new observational con ¬ 
stra ints for theoretical modelling. Recent work by Blana| 
et al.| ([2015]) raises questions about the current picture of 
the orbit of Hercules, and suggests that a path aligned 
along the major axis of the dwarf may not necessarily be 


correct. Our resul ts are consisten t with the orbital path 
proposed by Martin & Jin (2010), although both the po¬ 
sition angle of the major axis, and the alignment of the 
substructures suggest a slightly steeper orientation. In¬ 
terestingly, the position angle of the major axis that we 
determine does not line up perfectly with the alignment 
of the substructures (although the two are consistent 
within the 68% confidence band). This suggests a picture 
similar to that for the globular cluster Palomar 5, where 
the orbital path is no t quite aligned with the leading 
and trailing tidal arms (|Dehnen et al.|2004| |Odenkirchen 
et al.|[2009|. It may also be an indication that the tidal 


tails are transitioning to an orientation aligned towards 
the Mil ky Way as the dw arf approaches apo-galacticon 
(Klimentowski et al. 2009j). Finally, two of the significant 
over-densities we detect sit outside the confidence bands 
of the substructure alignment. This may suggest that 
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Hercules has made more than one peri-galactic pass, and 
is accompanied by free-floating debris. This would be 
consistent with infall timescales of 2-5 Gvr fou nd in the 
literature ( Rocha et al.|2012 Blana et al.|2015 ), and the 
destructiv e pen-galactic passage described by Martin fc] 


Jin 


( 2010 ). 

The identification of statistically significant over¬ 
densities provides new opportunities for spectroscopic 
follow up and further analysis in the ELT era, as well 
as new data that can be added to test and calibrate the¬ 
oretical models. 
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